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Abstract 



The exploration of atomic properties of strongly coupled partially degenerate plasmas, 
also referred to as warm dense matter, is important in astrophysics, since this thermodynamic 
regime is encountered for instance in Jovian planets' interior. One of the most important is- 
sues is the need for accurate equations of state and transport coefficients. The Ziman formula 
has been widely used for the computation of the static (DC) electrical resistivity. Usually, 
the calculations are based on the continuum wavcfunctions computed in the temperature- 
and density-dependent self-consistent potential of a fictive atom, representing the average 
ionization state of the plasma (average-atom model) . We present calculations of the electrical 
resistivity of a plasma based on the superconfiguration (SC) formalism. In this modeling, the 
contributions of all the electronic configurations are taken into account. It is possible to ob- 
tain all the situations between the two limiting cases: detailed configurations (a super-orbital 
is a single orbital) and detailed ions (all orbitals are gathered in the same super-orbital). The 
ingredients necessary for the calculation are computed in a self-consistent manner for each 
SC, using a density-functional description of the electrons. Electron exchange-correlation 
is handled in the local-density approximation. The momentum transfer cross-sections are 
calculated by using the phase shifts of the continuum electron wavefunctions computed, in 
the potential of each SC, by the Schrodinger equation with relativistic corrections (Pauli 
approximation). Comparisons with experimental data are also presented. 

1 Introduction 

The static electrical (DC) resistivity is important for the characterization of the plasma state. 
At high density, when the Spitzer theory becomes invalid, a consistent many-particle theory is 
required in order to evaluate the transport coefficients. The extended Ziman pQ formulation of 
the electrical resistivity of liquid metals is based on linear response theory, where free electrons 
in a metal are uniformly accelerated until they collide with an ion and are scattered. In the 
t— matrix formulation of Evans [2J, the resistivity reads [31 [HE]: 
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where = pAf/A is the matter density in cm -3 , p the matter density in g.cm -3 , A the 
atomic mass in g and TV the Avogadro number. The Fermi-Dirac distribution and its derivative 
read respectively: 

/(e,/x)= 1 + and ^feM) = -/3/(e,/x)[l-/(e,M)], (2) 

where (3 = 1/(/cbT) and p is the chemical potential. The function I is defined by an integral: 

1(e) = / S(q)a(q)q 3 dq, (3) 



o 



where k 2 = 2me/h 2 , m being the electron mass. The quantity S(q) is the ionic structure 
factor and a(q) represents the scattering cross-section, q 2 = 2k 2 (l — \) and x = cos(0), where 
6 = (k, k'), k and k' being the momenta of the conduction electron respectively before and after 
the scattering event. Therefore, one has: 



I(e)=2k 4 J ^S(V2k 2 (l- X ))\a(k,x)\ 2 (l-x)dx, (4) 

where a(k, x), defined by a(q) = \a(k, x)| 2 , is the scattering amplitude. In fact, the resistivity 
takes now a Drude-like form 

V = Vo 3— (5) 

n e % t 

with r/o = ha /e 2 =2l.74 /iO.cm, a =52. 91772083 10" 10 cm being the Bohr radius. n e = Z*rii 
represents the electron density. The characteristic time To = Aneohao/e 2 and the relaxation time 
r are such that 



7 = 3^ </>e WhCTe %(t^)mde. (6) 

Note that any model which treats the continuum by a Thomas-Fermi approximation will 
not give the correct continuum density of states and not correctly incorporate the effects of 
electron-electron interactions in the electron-ion phase-shifts [6j. Such a model is not reliable 
for the evaluation of the scattering cross-section. 



2 Electronic-structure model and electron-ion interaction 

The average-atom (AA) model we used relies on the following self-consistent calculation. Schrodinger's 
equation is solved in the Pauli approximation [8] , in which only first-order relativistic correc- 
tions are retained: 

Ve,l{r) = e y e ,ii r ), ( 7 ) 

where V mv is the mass-velocity term, Vd is the Darwin term, and V so is the spin-orbit 
contribution. y e ,r( r ) is the radial part of the wavefunction multiplied by r. The potential splits 
into two parts V sc {(r) = V e \(r) + V xc (r), where V xc represents the exchange-correlation part, 
evaluated in the local-density approximation [SJ, and V^i the electrostatic part. Outside the 



2m 



dr 2 



1(1 + 1) 



eV sci + VW + Vd + V sc 
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Wigner-Seitz (WS) sphere, i.e. beyond the radius r ws , a bound-electron wavefunction is given 
by: 



y n i(r) = A nl r Ki(-iKr), with K 



2me I < 



2Er 



1/2 



(8) 



where Ki is a modified spherical Bessel function of the third kind and A n \ a constant fixed 
by the boundary condition of the wavefunction at r ws and the normalization. Such a boundary 
condition will be noted BC2 in Sec. 01 while the usual "confining" boundary condition y n i(r ws ) = 
will be referred to as BC1. Beyond r ws , a free-electron radial wavefunction is written 



2mK ( e 
( 1 + E~o 



1/2 



r [cos[5i(k)] ji(Kr) - sm[5i(k)] ni(Kr)} , 



0) 



where ji and ni are the spherical Bessel functions of the first and second kinds respectively, 
Eq is the rest-mass energy of the electron, and 5i(k) is the phase shift, given by 



taaftCfc)] 



k j\ {kr v 



t,ji(kr 



(10) 



k n'jikryjs) - £ ni(kr ws ) ' 

where £ is the logarithmic derivative of the radial wavefunction y er (r)/r. The calculation of 
the phase-shifts Si(k) is painstaking and requires accurate numerical methods [TOj [HI E2] • One 
defines, for a fixed chosen A^k, the finite difference 



A^S k = max { Si(k) -Si(k- A^k) + 5 t (k + A^k) - S t (k) } . 

(i) 



(11) 



Then, if A^6 k > h (h being an infinitesimal threshold), iVf = AW^/Zj points are added 



in the interval \k — A^k, k + A^/c] , which enables one to define 

A^5 k = 5i(k) -Si(k- A^k) + 5i(k + A^k) -6i(k) 



(12) 

where A^k = 2A^k/N i k 1) . Then, for each k of the new grid, if A^5 k > h, then A^ 3) = 
A^5k/h points are added in the interval \k — A^k, k + A( 2 )/c] , etc. The iterative process goes 
on until the resonance is correctly described. The radial electron density is given by 



-2^+1)2/^) 



(13) 



z=o 



The new potential V^i, solution of Poisson's equation, is inserted in Eq. ([7]) and the process 
continues until convergence is reached. 

3 Scattering cross-section 

The electron-ion cross-section for momentum transfer reads 



a(q) = \a(k,xW 



1 



^(2i + l) e «'( fc )sin[5,(fc)]fl(x) 



1=0 



(14) 
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where P/(x) are the Legendre polynomials. It is useful to split the summation over I in 
two parts in Eq. ()14)) : the first summation running from to Iq (typically equal to 30) and 
the second one from (Iq+1) to oo. The latter can be approximated by the classical expression 
of the cross-section of a screened charge (summed from to oo) minus the values for I < Iq. 
Assuming, at high electron energies, a classical scattering cross-section for a screening potential 
a(q) ~ A 2 /(q 2 + A 2 ) 2 , this consists in writing [3] 



\ £ (2l + l)e i5 ^sm[S l (k)}P l ( X ) = + l)Qi(C)fl&) 



i=i +i 



2k 2 1 



Z=0 



lo 



-j2^+m(0Pi(x)) 



A 

2k 2 



C-x 



g(2l + l)Q,(C)fl(x) 



(15) 



with yl = 2mZ*e 2 /h 2 and C = 1 + X 2 /2k 2 , where A is a screening length. The quantity Qi(() 
represents Legendre function, evaluated using recursive relations for low x and from an integral 
representation (see Ref. [13], p. 1017, § 8.821) for x > 1.02. One gets 



a(k, X ) = lj2(2l + l)P l (x) 

1=0 



sin[<S,(*#* ,(fc) - Ag i(C) 



+ 



A 1 



Since e 2iS '^ 
the ratio 8i (k)/5, 



1 ~ 2i5i(k) = iAQi(z)/k [3], and since <^(fe) is known, for Z=0, • 
i -i(k) to Qi (z)/Qi _i(z) in order to determine z and then A 2 



(16) 



lo, one fits 



2A: 2 



1) 



In that way, A is a function of k. It can be numerically costly to solve this equation for the 
numerous values of k. In order to speed up the computation, one can use the asymptotic form 
Q n (z) » 2 n (n!) 2 /[(2n)!(2n + l)z n+1 ], which leads to Q n+1 (z)/Q n (z) w (n + l)/[(2n + 3)z], 
i.e., for n = 30, z ~ 0.492063 529/^30- The maximum energy of the continuum is taken to be 
£m = A* + 40 ksT. Therefore, at high temperatures, it can be relevant to add the following 
correction to the resistivity: 



At/ 



3ir e 2 n 2 r ' 



q 3 S(q)a(q)dq. 



(17) 



4 Comparisons with different theoretical models 

Rinker [14J proposed a combined structure factor which has the virtue to approach Debye-Hiickel 
limit as the core exclusion radius tends to 0, and the simple hard-sphere limit as the Debye radius 
tends to oo. The structure factor is expected to be important at low temperature only (T < 
10 eV according to [1]), but as can be seen on Table [H except for the Debye-Hiickel case, the 
differences between Percus-Yevick, BD (structure factor proposed by Bretonnet and Derouiche 
|15j for the OCP) and Rinker [T3] are quite small even at 5 eV. 

Table [2] shows the impact of the value of the bound wavefunction at r ws (fifth and sixth 
columns). The maximum difference can reach 25 %. In the case of iron at T=5 keV and /3=1965 
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p/po 


PY 


DH 


Rinker [14] 


BD [15J 


0.1 


732.1 


642.5 


731.6 


723.9 


0.2 


234.4 


158.0 


233.9 


224.3 


0.3 


212.8 


128.3 


212.0 


196.3 


0.4 


180.8 


101.1 


179.9 


163.4 


0.5 


142.6 


74.7 


141.8 


126.8 



Table 1: Resistivity r/ (/zfi.cm) of Al at T=5 eV for different structure factors: PY: Percus- 
Yevick (Hard Sphere), DH: Debye-Hiickel, Rinker (Ref. p3j, Eq. (15), p. 4211), and BD: 
Bretonnet-Derouiche [TS], structure factor for the one component plasma (OCP). 



p/po 


T(eV) 


Vps MS (HJ) 


r/p S AA 


r]BCi AA 


7]BC2 AA 


77 SC 






(01) 


(W) 


(This work) 


(This work) 


(This work) 


0.1 


5 


810 


405 


718.0 


731.6 


642.6 


0.2 


5 


720 


590 


231.9 


233.9 


316.1 


0.3 


5 


490 


318 


205.9 


212.0 


237.4 


0.4 


5 


312 


240 


132.0 


179.9 


200.2 


0.5 


5 


203 


165 


106.1 


141.8 


152.4 


0.2 


10 


382 


288 


296.1 


360.7 


375.9 


0.3 


10 


296 


240 


198.8 


200.5 


208.7 


0.4 


10 


220 


220 


151.8 


175.0 


185.0 


0.5 


10 


164 


164 


129.6 


149.2 


158.5 



Table 2: Resistivity r/ (^uQ.cm) of Al: comparisons with Table II of Perrot and Dharma-Wardana 
[3]. Their resistivities presented here use the pseudo-potential (ps) formulation. MS: Multiple 
Scattering and AA: Average Atom. BC1: the wavefunction cancels at r ws , BC2: the wave- 
function extends outside the WS sphere (jellium model, see Eq. ©). SC: Superconfiguration 
averaged. 



T=5 keV, /)=1965 g/cm 3 


Resistivity rj (/j,Q.cm) 


Rinker [14J A 


0.827 


Rinker [14J B 


0.649 


This work (BC1) 


0.693 


This work (BC2) 


0.701 



Table 3: Comparisons with theoretical results of Rinker [TJj for Fe. Model A corresponds to a 
Thomas- Fermi-Dirac potential and the combined structure factor of the third column of Table 
[TJ Model B uses a Thomas-Fermi-Dirac potential with a Debye-Hiickel structure factor. 
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Figure 1: Comparisons with calculations of Yuan et al. PS] and Perrot and Dharma-Wardana 
(Table I of Ref. [lj: pseudo-potential + scattering formulations). 
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Figure 2: Comparisons with the Benage experiment [19] for Al and with the calculations of 
Perrot and Dharma-Wardana [3], Lee and More [21] and Kitamura and Ichimaru [20J. 
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g/cm 3 (see Table ED, the difference between BC1 (the wavefunction cancels at r ws ) and BC2 
(which allows a bound-electron wavefunction to extend outside the WS sphere, see Eq. 
drops to 1.14 %. Our results lie between the results of the models A and B of Rinker [J4]. In 
the following, we keep the condition BC2 (atom in a jellium). 

We can see in Table [2] that our values decrease more rapidly with compression than the values 
of Ref. [1] (with or without multiple scattering), and are almost always smaller. The comparisons 
of Fig. Q] show that our results are closer to the pseudo-potential results including multiple 
scattering [3] than to the scattering- type calculations [H [16]. Some other approaches exist; 
for instance, Kuhlbrodt et al. [IT] calculated the resistivity within the linear response method 
of Zubarev. This correlation method is applicable for arbitrary degeneracy of the plasma, i.e 
from the low-temperature Spitzer regime up to the strongly coupled, degenerate domain where 
Ziman's approach is appropriate. It seems, however, that reliable resistivity calculations must 
include a unified description of the plasma composition and the density dependence of the inter- 
particle interactions |17[ [T5] . 

5 Comparisons with experimental values 

Resistivity measurements in plasma state are scarce so that the theoretical approach is in many 
cases the unique tool available in these studies. The measurement of resistivity along a thermo- 
dynamic path going from the melting point under normal conditions to a state corresponding 
to a compression 0.03 at T=26.5 eV in Al has been performed at Los Alamos by Benage et al. 
|19] using the exploded wire method. Figure [2] shows that a good agreement is obtained. The 
accuracy of the experimental results is generally within a factor of 2, which is of the same order 
of the typical discrepancy between theory and experiment. It is important to keep in mind that 
expression @ neglects multiple-scattering processes, which was shown to lead to a difference of 
a factor 1.5 in the theoretical values Our values are smaller than the results of Perrot and 
Dharma-Wardana [3] and the results of the density-response model of Kitamura and Ichimaru 
[20] . In the intermediate density range, our results agree fairly well with the analytical model 
of Lee and More [21] . which was rather unexpected. 

6 Effect of population fluctations 

The number of relevant bound-electron configurations in a plasma can be immense, especially as 
the atomic number Z increases. As a configuration is defined by orbitals occupied by an integer 
number of bound electrons, a superconfiguration (SC) is defined by super-orbitals, which are 
groups of orbitals close in energy (i.e. whose energies differ by less than a small fraction of fcsT) 
and that are populated in all possible ways consistent with the Pauli exclusion principle. For 
instance E = (Is2s2p) 10 (3s3p) 7 (3d) 2 (4s4p4d4/5s...) 1 is a SC composed by four super-orbitals. 
A configuration is a particular SC in which each super-orbital contains only one sole orbital. 
A reasonable number, i.e., a few hundred, of SCs can contain a large number of configurations 
(especially for high values of Z). The value of thermodynamic quantity X is given by 



where is the probability of the configuration 3 and Yls We = 1- A self-consistent 




(18) 
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calculation similar to the one described in Sec. [5] is performed for each SC. Table [5] displays the 
values of the resistivity averaged over the superconfigurations (third column) . 

A correction is required in the use of Z* instead of Z= since this occurs as Z% in the scattering 
cross-section. That is, we need to use (Z%) rather than (Zg) 2 . The fluctuations in Zg are 
intimately related to the fluctuations of the electron and ion densities of the plasma. Assuming 
that ion fluctuations are too slow to follow the electron-density fluctuations which determine 
the electronic configurations and electron-scattering processes, we consider the fluctuations in 
the electron subsystem only. We have 

< Zl >=< Z s > 2 (1 + nae/P), (19) 

where the finite-temperature interacting electron-gas compressibility Xe can be obtained from 
the electron-gas response function via the static structure factor since {Zg)riiXe = f3S ee (q = 0). 

Concerning the prefactor of the Ziman formula, since one has (l/Zi) > (1/Zg) 2 (Jensen's 
inequality for convex functions) and assuming that this dependence overcomes the one from the 
scattering cross-section, the resistivity averaged over SCs is expected to be larger than the one 
obtained from the AA model. However, it is important to keep in mind the fact that the latter 
prefactor should not be but rather (-^z)j^j, which could impact the results. This is 

related to Rinker's remark [22] that 1/Z* 2 is in fact \j(Z§Zi) where Zq can indeed be identified 
with Z* on variational grounds (Boltzmann equation), which is not the case of Zj, number of 
charge carriers, involved in the relaxation time. 

7 Conclusion 

We presented calculations of the DC electrical resistivity of warm dense plasmas using the 
generalized Ziman formula. The calculations do not use any model potential, dielectric function 
or structure data, but proceed via first-principle calculations based on density-functional theory 
in the local-density approximation. The electronic structure of an average scatterer is determined 
solving Schodinger's equation. The calculations require phase shifts which satisfy the finite- 
temperature Friedel sum rule. Different model structure factor have been tested, with two 
different boundary conditions of the wavefunctions. We discussed the approximations in the 
method arising from the use of an average-atom model and compared our results with other 
existing theoretical models. Our values are in fairly good agreement with experimental data. 
The method proposed here provides a self-contained and self-consistent approach to the electrical 
conductivity of plasmas of arbitrary degeneracy, structure and density. It enables one to perform 
calculations beyond the average-atom model, taking into account the fluctuations due to the 
electron configurations describing different ionization and excitation states. 
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